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Abstract — Main focus of this paper is control of infusion 
pump actuator (IPA) for biomedical applications. Control of 
IPA done using position regulator and control which combine 
position and speed regulator is proposed. The main guideline is 
to decrease power losses, heating and current consumption in 
the system. One of the important areas within the field of 
variable speed motor drives are the operational boundaries of 
system. Presently, the operational boundaries of variable speed 
motor drives are set based on the operational boundaries of 
single speed motors, i.e. by limiting current and power to rated 
values. This results in under-utilization of the system, and places 
the motor at risk of excessive power losses. In this paper, the 
first part represents general concept of infusion pump actuators 
and mathematical modeling of the system, then the control of 
position and speed of infusion pump actuator (IPA) by using 
fractional order PID controller is presented. The second part 
shows simulations of our own design of control speed and 
position by fractional order PID controller using SIMULINK 
and MATLAB software. Finally fractional order PID regulator 
results are compared with classical PID regulator results 

Index Terms — DC motor, infusion pump actuator (IPA), PID 
controller, FOPID controller, position control, speed control, 
combine position and speed regulator, simulations 


I. INTRODUCTION 

In reference to the literature ([1,2] and many others not cited 
here), for the better part of 20 th century most motion control 
systems were designed to operate at a fixed speed. Many 
existing systems still operate based on a speed determined by 
the frequency of the power grid. However, the most efficient 
operating speed for many applications, such as fans, blowers 
and centrifugal pumps, is different from that enforced by the 
grid frequency. Also, many high performance applications, 
such as robots, machine tools and the hybrid vehicle, require 
variable speed operation to begin with. As a result, a major 
transition from single speed systems to variable speed systems 
is in progress. 

The transition from single speed drives to variable speed 
drives has been in effect since the 1970s when movements 
towards conservation of energy and protection of the 
environment were initiated. About seventy percent of all 
electrical energy is converted into mechanical energy by 
motors in the industrialized world. This may be the most 
important factor behind today’s high demand for more 
efficient motion control systems. A large body of research is 
available on variable speed drives due to the significant 
industrial and commercial interest in such systems. 
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However, some areas of importance merit further 
investigation. One such area is the operational boundary of 
motor drives. The operational boundaries of variable speed 
motor drives are being incorrectly set based on the 
operational boundaries of single speed motor drives, i.e. by 
limiting current and power to rated values. The operational 
boundary of any motor drive must be set based on the 
maximum possible power loss vs. speed profile for the given 
motor. Also, all control strategies for a machine must be 
analyzed and compared based on such an operational 
boundary [3]. 

The mechatronic systems represent one of the most 
challenging control applications due to their interdisciplinary 
nature [4,5,6]. Numerous control algorithms have been 
proposed to deal with nonlinear Dynamics of the mechatronic 
systems. For linear mechatronic systems, the 
proportional-integral-derivative (PID) controller is often used 
owing to its simple structure and robustness [7]. Another 
approach in dealing with mechatronic systems challenges is 
the fractional-order (FO) control strategies. One of the most 
common applications in all mechatronic domains is the 
control of DC motors. The control of DC motors has been the 
interest of many researchers, due to the wide variety of 
applications that require the use of different types of DC 
motors [8, 9, 10]. The controllers designed for these DC 
motors range from simple traditional PIDs to advanced 
control algorithms, among which fractional order control has 
been gaining more and more popularity [11,12,13]. 

Fractional calculus has been used relatively recently in 
modeling and control applications [14,15]. The attractiveness 
of the fractional order PID controllers resides in their 
potential to increase the closed loop performance and 
robustness of the closed loop system, due to the extra tuning 
parameters available, as compared to the conventional 
controller. With fractional order controllers, the order of 
differentiation and integration may be used as supplementary 
tuning parameters and thus more specifications can be 
fulfilled at the same time, including the robustness to plant 
uncertainties, such as gain and time constant changes [16,17, 
18]. In general, frequency domain tuning of the fractional 
order controllers is preferred using optimization routines to 
yield the final solutions. The performance criteria are 
frequently specified in terms of gain crossover frequency, 
phase crossover frequency, phase margin, gain margin, 
robustness to open loop gain variations [19, 20]. 

However, some areas of importance merit further 
investigation. One such area is the operational boundary of 
motor drives. The operational boundaries of variable speed 
motor drives are being incorrectly set based on the 
operational boundaries of single speed motor drives, i.e. by 
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limiting current and power to rated values. The operational 
boundary of any motor drive must be set based on the 
maximum possible power loss vs. speed profile for the given 
motor. Another area is the regulator design for the speed and 
position. Also combine speed and position regulator design is 
to impose some additional control on the same system. 
Recently researches for speed and position regulators are 
focused on fractional order PID and classical PID regulators 
[21,22]. Therefore combine speed and position control of 
IPA is design and tested by using fractional order PID as well 
as classical PID 

The paper is structured as follows: Section II introduces 
general concept and IPA mathematical model, section III 
presents the tuning methods for the fractional order PID 
controller, with the position control for a DC motor. Section 
IV presents MATLAB design of combine speed and position 
regulators with fractional order PID and classical PID. 
Section V presents SIMULINK and MATLAB simulations, 
while the last section includes the final conclusions. 

II. GENERAL CONCEPT AND MATHEMATICAL 
MODEL OF IPA 

Linear actuator is a device that applies force in a linear 
manner, as opposed to rotationally like an electric motor. 
Electro-mechanical actuators are similar to mechanical 
actuators except that the control knob or handle is 
replaced with an electric motor. Rotary motion of the 
motor is converted to linear displacement of the actuator. 

A. Mathematical Model and Control of Inf ision Pump 
System 

a) DC Motor Speed and Position Control 
Ra La > 

+ 


Fig 1. DC motor model. 



A general model of a DC motor is shown in Figure 1. Applied 
voltage V a , which is the manipulated variable, For the 
position conrol 0(t), which is the control variable. For the 
speed control, the controlled variable is the angular velocity 
w(t) and the load transfer function has the form in (1): 


G(s) 


W(b) 

Va(e) 


Km 

(Las+RsXjB-t-bJ+KbKm 


( 1 ) 


However, for many dc motors the time constant of the 
armature x a =L a /R a is negligible and therefore the model can be 
simplified to : 


G^} = 


is+i 


( 2 ) 


goo = m = _^ 

Va [tfi a (Ta+ 1) 


( 3 ) 


Injection mg is calculated based on 10 cc injector. The dose 
formula for the infusion pump mechanism shown in Figure 2 
is given by (4). 

Dose (mg) = (S / 2n) a.m'} n d ilaf (4) 


Where, the angular position of the 0 (rad) shows the screw 
density, the screw pitch a (mm), the injector radius r en _ (mm) 
and the drug concentration d (mg / mm 3 ). 



Fig 2. Syringe infusion mechanism [8]. 

In this study, a syringe mechanism was used to serve as the 
fluid reservoir and to create the necessary pressure for 
infusion. It uses an endless worm thread to position the 
plunger, managing in this way the liquid movement. This 
process is controlled by a dc motor where the rotation of the 
motor is transmitted to the endless thread (Figure 1. ). 
Normally, a spring is used to push the plunger with a constant 
force, necessary for a stable infusion pressure. This type of 
infusion generates a continuous high precision flow (with an 
error below 2%). 

III. FRACTIONAL ORDER PID CONTROLLER DESIGN 


The fractional order PI^D^ controller [15, 16] was proposed 
as a generalization of the PID controller [2] with integrator of 
real order X and differentiator of real order (i. The transfer 
function of such a type controller in the Laplace domain takes 
the form 

G c (s)=K p + K lS -* + K d ^ (5) 

Here k p , ki , k d denote the proportional, integral and 
differential gains of the controller respectively. k p fi and k d 
parameters usually take non-negative values. Taking k = 1 
and /u=l we obtain a classical PID controller. If k d = 0 or k t = 
0 we obtain PI A controller or PD^ controller, respectively. All 
these types of controllers are particular cases of the PI^D^ 
controller. In this paper we use the both PD^ and PI A 
controllers, whose transfer function have the form, 
respectively 

i j c (s) = -F kg , G c (s) = k p -F k { s~ l (6) 


Where x=JR a /(R a b+K b K m ) and K dc =K m / (R a b+K b K m ) 

The transfer function from position O(t) as output (controlled 
variable) to armature voltage V a as input (manipulated 
variable) will be : 


The block diagram of the considered feedback control system 
is shown in Figure 3, where G p (s) has form (1) and G c ( l s') has 
form (5). 
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Fig. 3. Feedback control system structure 

A. FOPID Controller Design in Frequency Domain 

In the frequency domain design method for fractional areas, 
the controllers first start with the linear model of the system. 
The transfer function of the G p (s) and G c (s) controllers 
supports the properties of the following equations. 

a) An imposed phase margine of the open loop system 


parameters can be specified. The stability boundaries are 
curves on which each point corresponds to polynomial (4) 
having zeros on the imaginary axis. It may be the real zero 
boundary or the complex zero boundary. It is easy to see that 
polynomial (4) has zero s = 0 if k p = 0 (the real zero 
boundary). The complex zero boundary corresponds t8o the 
pure imaginary zeros of (12). We obtain this boundary by 
solving the equation 

w(jco) = e~ J *K(k p + k d (ja>Y)+ jamjar + 1) = 0, (13) 

which we get by substituting s = join polynomial (12) and 
equating to 0. The complex Equation (13) can be rewritten as 
a set of real equations in the form 


arctan ( G.j, (jtd’) G p (jtiO J = —it -f pm (7) 

b) An imposed gain crossover frequency of the open 

loop system 

|G c (fu)G p (jw}| - 0 dB (8) 

c) A condition for robustness to gain variation 

arctan ( g c Cjkis) G p (jkia) J = 0 (9) 

d) An imposed remove high frequency noise 

U>Ll4iXl U n l.ltijJ 

= B dB ( 10 ) 

e) An imposed remove the output distortion 


AaB (l 1 

Controller saturation determines the physical systems of the 
agents. For the motor-generator system, the control operation 
is limited to 10 V. All equations should be used from the 
equation (7) to the equation (10) to design features such as 
phase margin, frequency shift gain, amplitude of sensitivity 
function (A) and size of complementary sensitivity function 
(B) for efficient desired design features of the FOPID 
controller 


B. Parametric Synthesis of FOPID Controller 

The asymptotic stability is the basic requirement of a 
closed-loop system. Gain and phase margins are measures of 
relative stability for a feedback system, therefore the synthesis 
of control systems is very often based on them. In typical 
control systems the phase margin is from 30° to 60° whereas 
the gain margin is from 5 dB to 10 dB. 

In this paper we tuning the controller for the specified phase 
margin, because a phase margin is closely related to a 
transient response, i.e. overshoot. For this purpose we use the 
phase margin tester exp(~ 7 ^), where ^is a phase margin. This 
tester does not exist in the real control system, it is only used 
for the controller tuning. 

Taking into account the phase margin tester exp(-y^) in the 
main path of control, plant (l)and controller (3) we obtain the 
characteristic polynomial of the closed-loop system 


Re[ w(jco)] = 0, Im[ w(jco)] =0. (14) 

where Re[w(y co)] and Im[w(y co)] denote the real and the 
imaginary parts of (13), respectively. Finally, by solving the 
Equations (14) we obtain 


K 


1 _ rt wwin.0-™s0 

= U - 1 P . ,.W. 

JfSLfi j -f) 


vei sin! 0) 


(15) 

(16) 


Equations (15) and (16) determine the complex zero 
boundary as a function of co. The real zero boundary and the 
complex zero boundary for co > 0 decompose plane (k h k p ) 
into some regions. The stability region is chosen by testing an 
arbitrary point from each region and checking the stability of 
polynomial (12). In this paper only the stability region in the 
parameter plane of polynomial (12) is presented. When the 
stability regions are known, the tuning of the fractional 
controller can be carried out. 


Taking into account the phase margin tester exp(-y^) in the 
main path of control, plant (l)and controller (3) for speed 
control we obtain the characteristic polynomial of the 
closed-loop system 

w(s) = e~ 1 *K(k p + k i s^) + ST +1 (17) 

Using the classical D-partition method the stability region in 
the parameter plane (k i9 k p ) may be determined and the 
parameters can be specified. The stability boundaries are 
curves on which each point corresponds to polynomial (17) 
having zeros on the imaginary axis. It may be the real zero 
boundary or the complex zero boundary. It is easy to see that 
polynomial (4) has zero s = 0 if k p = 0 (the real zero 
boundary). The complex zero boundary corresponds to the 
pure imaginary zeros of (17). We obtain this boundary by 
solving the equation 

w(jco) = e~ j</> K(k p + k t ( jco ) ~ x ) + jcor +1= 0, 

(18) 


w(s) = e~ j</> K(k p + k d s M ) + s(st +1). (12) which we get by substituting s =jco in polynomial (9) and 

equating to 0. The complex Equation (10) can be rewritten as 
Using the classical D-partition method the stability region in a set of real equations in the form 
the parameter plane (k d , k p ) may be determined and the 
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Equations (12) and (13) determine the complex zero 
boundary as a function of co. The real zero boundary and the 
complex zero boundary for co > 0 decompose plane (k h k p ) 
into some regions. The stability region is chosen by testing an 
arbitrary point from each region and checking the stability of 
polynomial (9). In this paper only the stability region in the 
parameter plane of polynomial (9) is presented. When the 
stability regions are known, the tuning of the fractional 
controller can be carried out. 

A MATLAB GUI programs are developed the selection of 
control parameters for FOPID and PID controllers. Program 
outputs are output screen so that users enter the control 
parameters of PID controller as well as FOPID controllers. 
Control parameter values observed on the screen is obtained 
for dc motor speed control transfer function (l).Also this 
programs are used for dc motor position transfer function (2) 
by tuning control parameters of FOPID or PID controller. 


Fig. 5 MATLAB Block diagram of DC motor and combined 
position and speed regulator 

V. COMPUTER SIMULATIONS 

In this section, SIMULINK is used in the design of fractional 
order PID controller for both position and velocity controller 
of IPA and for modeling of IPA system. In order to compare 
the performance of classical PID and fractional order 
controllers two types SIMULINK diagrams are presented 
separately for both classical and fractional order controllers. 
In general PD controller is used for position control systems 
and PI controller is used for velocity control systems. In 
simulations both classical and fractional order PI and PD 
controllers are implemented 

A. FOPID Simulink Diagram 


IV. COMBINE POSITION AND SPEED PID 
REGULATOR 

Main guideline is to impose some additional control on the 
same system. Position PID regulator will stay unchanged. 
Imposed component is output signal from the speed controller 
shown in Figure 4. Speed controller will take effect on the 
whole system like additional differential component so 
coefficients in speed controller need to be the same order like 
differential component in position controller. Combined 
regulator shown in Figure 5 will be based on sum of responses 
of both position and speed regulators. Therefore combine 
regulator use two type PID controllers. PI controller is 
implemented for speed control and PD controller is used for 
position control. Also we can design combine position and 
speed fractional order PID regulator by implanting two 
fractional order PID controllers. 



Fig 4 Block diagram of DC motor and position(speed) PID 
regulator 


FOPID control parameters are determined from the speed 



transfer function as on screen shown in Figure 6 . Also FOPID 
parameters can be similarly find for the position transfer 
function by using same method. Screen values for speed 
control is depicted in Figure.6 


Fig 6 . FOPID control parameter values for speed control 
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FOPID SIMULINK diagram for FOPID position and speed 
regulators is shown in Figure 7. In this diagram there are one 
input (reference input is position) and two outputs, position 
and velocity respectively. Velocity and position waveforms as 
a simulation results are depicted in Figure 8a, 8b.respectively. 




Fig 9 . Classical PID control parameter values seen from 
output screen for speed control 


Fig 7. Fractional order PID SIMULINK diagram 



Time [s] 

Fig 8a. FOPID velocity waveform of dc motor 



Fig 8b. FOPID Position Waveform Of Dc Motor 
B. PID Simulink Diagram 

MATLAB -GUI screen shown in Figure 9, which is 
developed in section 4 is used to choose control parameters 
for speed transfer function. PID control parameter values on 
the screen is choosen for GOO =-=^7 speed transfer 

1.45.3 +i r 

function 


PID SIMULINK diagram for PID controllers is shown in 
Figure 10. In this diagram there are one input (reference is 
position) and two outputs, position and velocity respectively. 
Velocity and position waveforms as a simulation results are 
depicted in Figure 11a, lib respectively. 



Fig 10. Classical PID SIMULINK diagram. 



Tine [s] 

Fig 11a. Classical PID velocity waveform of dc motor. 
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Time [s] 


Fig lib. Classical PID position waveform of dc motor 

The Simulation results revealed good performances and show 
a stable and convergent behavior of the system when dealing 
with fractional order control law. The performances of both 
fractional order PI and PD controllers are analyzed and 
compared with integer order PI and PD controllers. The 
simulation results show that the fractional order controllers 
outperform the classical integer order controllers. 

VII: CONCLUSION 

In this paper, a design of a fractional order PI and PD 
controllers has been made to control the speed and position of 
a DC motor. The fractional order controllers were designed 
based on frequency domain specifications. The Simulation 
results revealed good performances and show a stable and 
convergent behavior of the system when dealing with 
fractional order control law. The performances of both 
fractional order PI and PD controllers are analyzed and 
compared with integer order PI and PD controllers. The 
simulation results show that the fractional order controllers 
outperform the classical integer order controllers. 

Main guideline is to impose the some additional control on the 
same system. Position PID regulator will stay unchanged. 
Imposed component is output signal from the speed 
controller. Speed controller will take effect on the whole 
system like additional differential component so coefficients 
in speed controller need to be the same order like differential 
component in position controller. Combined regulator will be 
based on sum of responses of both position and speed 
regulators. 

Having in mind obtained results, one can conclude that: 

• Small signal bandwidth of both types of regulation will be 
the same but response of combined regulator will be 
smoother; 

• On maximum requirements, combined regulation can 
produce more overshoot, but in biomedical application that 
type of signal will never be used in real situation; 

• Combined regulation generates less noise for lower level of 
errors; 


• Results show that by using combined controller we can 
reduce current consumption from battery. That will have as a 
benefit less heating in the system and possibility of longer 
flight for the same power supply in one system; 

• Control alternation (change of type of control) could be done 
in software code - no hardware change needed. 
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